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DOPPLER  PROPAGATION  STUDY  pages 

31  January  1963  Contract  AF04(695)-113 

This  report  describes  Doppler  propagation  and  the  problems 
encountered  in  transmitting  and  receiving  Doppler  tracking 
signals  through  the  atmosphere.  It  presents  formulas  for 
the  Doppler  shift,  models  of  the  atmosphere,  and  equations 
for  atmospheric  refraction.  The  information  is  combined  in 
an  analysis  of  Doppler  errors  caused  by  propagation  effects 
and  indicates  the  accuracies  to  be  expected  when  dealing 
with  Doppler  shifts. 


This  unclassificd  abstract  is  pssibnep  for  retention  in  a 
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LINE.  (IF  THE  ABSTRACT  IS  CLASSIFIEC,  HOWEVER,  IT  MUST  NOT 
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FOREWORD 


this  XeGhniGal  DoGumentary  Report  on  Definitive  Contraet  AF04(695)- 
113  has  been  prepared  in  aGGordanoe  with  Exhibit  A  of  that  contract  and 
Paragraph  3.18  of  AFBM  Exhibit  58-1,  "Contractor  Reports  Exhibit," 
dated  1  October  1959,  as  revised  and  amended. 

this  report  was  prepared  by  Philco  Western  Development  Laboratories 
in  fulfilling  the  requirements  of  Part  XVIII,  Schedule  of  Contract 
AF04(695)-113. 
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SECTION  1 
INTRODUCTION 

In  recent  times,  the  Doppler  effect  has  been  used  to  measure  the 
velocity  of  an  aircraft  or  satellite  with  respect  to  an  observer  or 
tracking  complex.  The  purpose  of  a  tracking  complex  is  to  character¬ 
ize  the  trajectory  of  a  vehicle  to  such  an  accuracy  that  reliable  pre¬ 
dictions  of  its  past  and  future  positions  can  be  made.  The  reliability 
of  a  particular  prediction  depends  upon  the  following  factors: 

a.  The  uncertainty  associated  with  the  satellite's  position  and 
velocity  vector  during  the  time  period  immediately  preceding 
prediction. 

b.  The  accuracy  of  the  mathematical  model  representing  the  satel¬ 
lite's  physical  environment. 

It  is  the  purpose  of  this  report  to  discuss  these  problems  with  respect 
to  the  Doppler  system.  The  Doppler  system  itself  is  described,  and 
equations  are  developed  for  the  Doppler  shift.  The  uncertainty  in 
vehicle  position  is  developed  following  a  discussion  of  the  mathematical 
model  representing  the  satellite's  physical  environment  (atmospheric 
characteristics  and  related  refraction  effects) . 
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SECTION  2 
DOPPLER  EFFECTS 

2.1  GENERAL  DOPPLER  EFFECTS 

ie 

The  Doppler  effect  is  observed  when  a  source  of  radiation  is  in 
motion  with  respect  to  an  observer.  It  is  characterized  by  an  apparent 
change  in  the  frequency  of  the  emitted  radiation.  Assuming  that  the 
observer  either  knows  or  has  a  method  of  monitoring  the  original  tranS* 
mitted  frequency  of  the  radiation,  the  relative  motion  of  the  radiating 
source  with  respect  to  the  observer  can  be  determined.  When  the  radi¬ 
ating  source  and  the  observer  are  moving  toward  each  other,  the  observer 
notes  a  "shift"  to  a  higher  frequency  caused  by  the  "foreshortening"  of 
the  radiated  wavelength.  When  the  source  and  observer  are  moving  away 
from  each  other,  a  decrease  in  frequency  is  noted,  caused  by  the 
"lengthening"  of  the  radiated  wavelength.  The  amount  of  frequency  or 
wavelength  shift  corresponds  to  the  relative  velocity  of  the  source, 
whereas  the  increase  or  decrease  in  observed  frequency  or  wavelength 
corresponds  to  the  relative  direction  of  movement  of  the  source. 

Doppler  data  systems,  for  example,  extract  frequency  shift  inform¬ 
ation  and  resolve  it  into  relative  velocity  measurements.  An  angle 
generating  system  is  frequently  coupled  with  a  Doppler  data  system  to 
provide  coordinated  angles  and  range  rates  for  the  purposes  of  orbital 
operations  and  ephemeric  constructions. 

2.2  ONE-WAY  DOPPLER  SHIFT 

Suppose  a  satellite  moving  with  a  velocity  of  magnitude  V  carries 
a  radio  transmitter  operating  at  a  frequency  f^.  The  Doppler  frequency 
shift,  f^,  is  the  difference  between  transmitted  frequency,  f^,  and 
received  frequency,  f^,  at  an  observation  point: 


*  Named  for  C.  J.  Doppler,  physicist  and  mathematician,  who  first  used 
this  effect  in  1842  in  determining  relative  velocities  of  stars, 
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this  Doppler  shift  may  also  be  measured  as  the  time  rate  of  change 
of  the  phase  difference  between  transmitted  and  received  signals,  the 
phase  difference,  0,  is  represented  by: 

t  .  s/c 

where  g  is  the  speed  of  light  in  a  vacuum  and  S  is  the  geometric  path 
length  between  transmitter  and  receiver,  the  optical  path  length,  S', 
may  be  expressed  as  an  Integral  involving  the  refractive  index: 

S'  •=  /  n(s)ds. 

0 

it  has  been  found  (Ref.  6)  that  the  following  inequalities  hold  between 
the  distances  D,  S,  and  S': 

D  <  S  <  S' 

S  -  D  <  S'  -  S 

where  D  is  defined  as  the  geometric  length  of  a  straight  line  from 
transmitter  to  receiver. 

It  is  important  to  recognize  that  the  refraction  error  in  range, 

6  R,  is  S'  -D,  not  S  -  D.  The  difference  between  S  and  D  is  usually 
negligible,  as  can  be  seen  by  approximating  the  curved  path  between 
transmitter  and  receiver  by  a  circular  arc  spanning  chord  D  and  sub¬ 
tending  a  central  angle  of  ££•  The  angle  CE  is  the  "refractive  error" 
in  the  apparent  elevation  angle  of  the  satellite.  This  yields  the 
result 
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or  10  ^ 
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for  most  elevation  angles  of  practical  use,  or  a  maximum  range  error  on 
the  order  of  feet  at  ranges  within  which  significant  refraction  errors 
will  occur.  The  phase  difference  then  becomes 

f  /? 

Q  s  ^  /  n(s)ds.  (2^1) 

'^0 

TO  get  the  time  rate  of  change  of  the  phase  difference,  the  above 
expression  is  operated  on  by  the  Leibnitz  rule  for  differentiating  an 
integral  involving  a  parameter: 


=  f^  (Doppler  Shift)  (2-2) 

where  n^  is  the  local  index  of  refraction  at  the  transmitter  correspond¬ 
ing  to  the  frequency  f^.  The  first  term  may  be  expanded  to  show  the  re¬ 
sults  of  three  effects.  The  first  is  the  actual  Doppler  frequency.  The 
second  is  an  error  term  due  to  refractive  bending  of  the  wave  path.  Such 
bending  results  in  a  direction  of  arrival  of  the  signal  at  the  vehicle 
different  from  the  straight  line  path  from  the  ground  to  the  vehicle. 

This  causes  the  measured  Doppler  velocity  to  differ  from  the  actual 
Doppler  velocity.  The  third  effect  in  the  expanded  first  term  in 
equation  (2-2)  is  an  error  arising  from  the  fact  that  the  phase  velocity 
of  the  wave  at  the  vehicle  is  not  equal  to  the  speed  of  light. 

The  last  term  in  equation  (2-2)  is  a  combination  of  two  types  of 
errors.  One  part  of  the  term  takes  into  account  the  error  due  to  a 
change  in  the  propagation  path  as  the  vehicle  moves  through  different 
strata  at  a  constant  radio  distance  from  the  ground  station-  The  other 
error  arises  from  irregularities  in  the  layers,  such  as  rain,  clouds, 
ion  concentrations,  etc.  Usually  the  second  term  in  equation  (2-2)  is 
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not  considered  or  measured  separately,  with  the  result  that  it  appears 
as  a  "noise"  or  "scintillation"  effect  riding  on  the  first  term.  Gen¬ 
erally,  it  is  neglected*  and  equation  (2-2)  reduces  to 

f  . 

t  Q  hs 

d  ”  c  \  dt  ' 


A  general  formula  for  the  Doppler  shift  in  one-way  transmission 
may  therefore  be  written: 


f 


d 


f 

o 


f 

r 


o  ds 
c  ^'t  dt 


or  alternately, 

‘oG  (2-3) 

A  sample  graph  of  one-way  Doppler  shift,  slant  range,  and  elevation 
angle  as  a  function  of  time  for  a  typical  satellite  pass  is  presented  in 
Fig.  2-1.  The  calculation  for  this  graph  assumed  a  spherical  earth,  a 
transmitted  frequency  of  400  Me,  and  a  vehicle  altitude  of  160  n.m. 

2.3  TWO-WAY  DOPPLER  SHIFT 

Consider  a  case  where  the  signal  is  transmitted  from  the  ground 
station,  received  in  Che  satellite,  and  retransmitted  to  the  ground 
station.  If  f^  is  the  frequency  transmitted  from  the  station  and  n^ 
is  the  index  of  refraction  at  the  station  corresponding  to  f^>  then 
the  frequency  received  at  the  satellite  is 

-  -  1  da  \ 

r  o  \  c  8  dt  / 


*  This  is  Justified  statistically.  A  discussion  of  the  subject  is 
given  in  Ref.  11. 
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the  ttransforined  frequency  in  the  satellite  Is 


f.  *a+bf  =  a+bf  ll-^n  4r 
t  r  o  \  G  8  dt 


(2^4) 


where  a  and  b  are  known  constants  which  depend  upon  the  equipment  in  the 
satellite.  (For  example,  a  =  0  and  b  «  1  in  the  special  case  of  skin 
tracking.)  If  n^  is  the  index  of  refraction  at  the  satellite  correspond¬ 
ing  to  f  ,  then  the  frequency  f  received  at  the  ground  Station  is 


g  t 


1  1  ds 

c  ’^t  dt 


(2-5) 


The  composite  Doppler  shift  is  found  by  combining  equations  (2-4)  and 
(2-5): 


f.  =  f  -  f  =  f  -f  (l-  —  n 
d  o  g  o  t\  ctdty 

<  r  1  SI  1 

«f  -  a+bf  11-  —  n  ,11 - 

o  o\  csdtyV  ctdty 

=  f^  (1  -  b)  -  a  +  -i  ^  [(a  +  byn,.  +  bf^nj 

Kf  Y 

o^s^t  \c  dt  J  ' 

( 1  dsV 

\c  dty  * 


Neglecting  the  last  term  because  of  the  smallness  of 
comes 


-  to  (1  .  b)  -  a  +  i  If  [(a  +  bfo)a^  +  btoHj  . 


this  be-' 


(2-6) 


In  the  case  of  skin  tracking  (a  «  0,  b 

-  ,  1  ds  ,  ^  V 

f  j  »  f  “  -JT  (n  +  n..) . 
d  o  c  dt  8  t 


1), 


(2-7) 
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As  previously  mentioned,  the  general  formulas  which  have  been  used 
here  are  affected  by  refractive  bending.  In  order  to  determine  the 
effect  of  this  on  accuracy,  the  next  several  sections  will  deal  with  the 
refraction  errors  in  elevation  angle  and  the  mathematical  model  of  the 
satellite's  physical  environment,  the  atmosphere.  The  equation  for  re¬ 
fractive  bending  will  be  derived  and  applied  to  correct  the  Doppler 
measurements. 
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SECTION  3 

ATMOSPIffiRlC  CHARACTERISTICS 


3 . 1  GENERAL 

As  a  radio  wave  passes  between  the  satellite  and  the  ground,  it 
is  affected  by  the  medium  through  which  it  travels.  The  propagation 
medium  will  attenuate,  reflect,  and  refract  the  wave.  At  the  frequen* 
Cies  under  consideration  for  tracking,  attenuation  and  reflection  by 
the  propagating  medium  are  of  a  negligible  magnitude;  but  the  effect 
of  refraction  normally  may  not  be  ignored. 

Atmospheric  refraction  deteriorates  the  accuracy  of  Doppler  measure¬ 
ments.  To  some  extent,  the  refractive  effects  can  be  predicted  and  the 
tracking  data  can  be  corrected.  As  the  various  portions  of  the  wave- 
front  encounter  atmospheric  regions  of  differing  refractive  indexes, 
their  velocities  will  be  increased  or  decreased  so  that  the  net  effect 
is  a  bending  of  the  ray  path  toward  the  perpendicular  in  regions  of 
higher  index  and  away  from  the  perpendicular  in  regions  of  lower  index. 
Various  methods  for  compensating  refractive  errors  have  been  devised, 
but  it  is  obvious  that  a  complete  correction  requires  a  continuous  pro¬ 
file  of  refractive  index  along  the  signal  path. 

3.2  INDEX  OF  REFRACTION 

The  index  of  refraction,  n,  for  electromagnetic  radiation  at  a 
given  point  in  the  atmosphere  is  defined  as  the  ratio  of  the  speed  of 
radiation  in  a  vacuum  to  the  speed  of  radiation  at  the  given  point. 

Since  the  atmospheric  index  may  differ  from  unity  by  several  hundred 
parts  per  million,  another  parameter  called  the  "refractivity"  is  often 
used.  This  parameter,  N,  is  numerically  more  convenient  and  is  defined 
as  follows: 


N  =  (n  -  1)  10^. 
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3.3  SPHERICALLY  STRATIFIED  ATMOSPHERE 

As  a  first  approximation,  the  shape  of  the  earth  is  considered 
spherical.  In  the  atmosphere  surrounding  this  sphere,  the  mean  density 
of  the  gas  decreases  approximately  exponentially  with  increasing  altitude. 
Usually  under  relatively  quiet  atmospheric  conditions,  the  vertical 
gradient  of  refractivity  is  much  greater  than  the  horizontal  gradient. 

(In  fact,  it  is  felt  that  the  azimuthal  or  horizontal  errors  in  Ref.  11 
are  less  than  10  percent  of  the  magnitude  of  the  elevation  or  vertical 
errors.)  For  Such  reasons,  but  principally  for  mathematical  simplicity, 
it  is  Customary  to  assume  that  the  refractivity  is  solely  a  function  of 
altitude  or  of  distance  from  the  center  of  the  earth.  However,  it  is 
understood  that  such  a  spherical  model  deviates  from  the  actual  condi¬ 
tions  in  at  least  two  ways.  Because  of  diurnal,  meteorological,  and 
seasonal  fluctuations  of  atmospheric  conditions,  the  refractivity  pro¬ 
file  varies  considerably  with  geographic  location.  In  some  cases,  the 
solar  zenith  angle  (i.e.,  local  time  and  geographic  position)  also  be* 
comes  an  important  factor  (Appendix  A  of  Ref.  3).  Further,  the  figure 
of  the  earth  is  ellipsoidal  rather  than  spherical. 

The  principal  differences  between  refraction  in  the  concentric 
spherical  atmospheric  stratification  and  in  the  confocal  ellipsoidal 
stratification  may  be  noted  as  follows.  The  angular  refraction  in  the 
zenith  direction  is  zero  in  the  spherical  model,  but  generally  non-zero 
in  the  ellipsoidal  case.  In  a  spherical  model,  the  ray  path  is  a  plane 
curve  and  its  angle  of  arrival  at  the  station  is  independent  of  the 
azimuth  of  the  ground  projection  of  the  ray.  In  the  ellipsoidal 
stratification,  the  ray  path  is  a  space  (nonplanar)  curve  and  its  angle 
of  arrival  at  the  station  depends  upon  the  azimuth  of  the  ground  pro* 
jection  of  the  ray.  For  instance,  considering  rays  which  make  equal 
45°  angles  with  a  normal  to  the  earth- spheroid  at  infinity,  the  ones 
reaching  the  station  from  the  east  or  west  will  arrive  at  an  elevation 
angle  slightly  greater  (about  2  seconds  of  arc  at  a  station  on  the 
equator)  than  that  of  such  rays  arriving  from  north  or  south.  Consid¬ 
ering  rays  which  are  inclined  at  45°  to  the  polar  and  to  the  equatorial 
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radii  at  infinity,  the  elevation  angle  of  arrival  at  the  equator  may  be 
some  11  minutes  of  arc  greater  than  the  corresponding  angle  at  the  pole 
due  to  the  difference  in  earth- curvature  alone.  Such  conclusions  de¬ 
pend  upon  ellipsoidal  geometry  and  the  nature  of  the  equations  for  re¬ 
fraction  in  the  concentric  spherical  stratification.  By  differentiating 
these  equations,  it  may  be  shown  that  the  local  mean  radius  of  earth- 
curvature  need  be  known  only  with  an  accuracy  comparable  to  the  inherent 
angular  accuracy  of  the  cracking  system  when  approximating  the  ellipsoidal 
stratification  by  a  spherical  stratification.  Since  the  inherent  angular 
accuracy  is  on  the  order  of  one  part  in  10^  in  the  most  accurate  present 
systems,  this  figure  represents  the  approximate  maximum  accuracy  for 
specifying  the  mean  radius  of  earth  curvature  at  a  tracking  station. 

Because  of  the  small  eccentricity  (0.082)  of  the  earth  spheroid, 
the  maximum  difference  in  the  principal  radius  of  curvature  at  the  Sta¬ 
tion  is  about  one  part  in  10^.  Hence,  in  making  refraction  corrections, 
the  model  of  a  spherically-stratified  atmosphere  is  sufficient  for  all 

existing  tracking  systems,  provided  the  local  mean  radius  of  earth 

$* 

curvature  is  determined  with  an  accuracy  of  one  part  in  10  .  The  cur¬ 

rent  accuracy  in  the  radius  is  close  to  this  value. 

In  considering  the  spherical  stratification  of  the  earth's  atmos¬ 
phere,  certain  layers  distinguish  themselves  one  from  the  other,  especial¬ 
ly  with  regard  to  index  of  refraction.  For  example,  the  lower  level 
(troposphere)  has  a  positive  refractivity  and  the  upper  level  (ionosphere) 
has  a  negative  refractivity. 

*  In  computing  the  mean  radius  at  a  station,  the  geodetic  latitude 
must  be  known  to  4  5  minutes  of  arc.  The  mean  radius  is 

e  b^(l  -  e^  sin^ji) 

where  ^  «  geodetic  latitude  of  the  station, 

e  s  eccentricity  of  the  earth  model,  and 
b  s  semi -minor  axis  (polar  radius) . 
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3.4  TROPOSPHERE 

The  troposphere  is  that  portion  of  the  atmosphere  below  about  40 
nautical  miles  in  which  clouds  form,  convective  disturbances  occur,  and 
temperature  decreases  with  altitude.  In  the  troposphere,  the  refrac- 
tivity  is  a  function  of  meteorological  variables  Such  as  water  vapor 
pressure,  air  temperature,  and  air  pressure  although  independent  of 
radio  frequency  up  to  30  kMc.  The  refract ivity  is  approximated  by 

N  =  (n  -  1)  10^  -  I  (P  +  1^) 

where  N  *  refractivity , 

n  =  index  of  refraction, 

T  =  temperature  in  degrees  Kelvin, 

P  =  air  pressure  in  millibars, 
e  s  water  vapor  pressure  in  millibars, 
a  =  constant,  77.6°  K/mb,  and 
b  =  constant,  4810°  K. 

Since  it  is  not  usually  convenient  to  determine  the  pressure  and 
temperature  at  an  arbitrary  point  in  space,  models  are  employed  to  rep¬ 
resent  the  troposphere.  Experimental  evidence  has  substantiated  that 
in  the  troposphere  the  mean  refractive  index  decreases  approximately 
exponentially  with  increasing  altitude,  while  also  varying  with  geo¬ 
graphic  location.  The  CRPL  Exponential  Radio  Refractivity  Atmosphere 
or  "exponential  reference  atmosphere"  is  a  good  single- function  rep¬ 
resentation  of  refractivity  for  the  troposphere.  The  model  calculates 
the  refractivity  as  a  function  of  general  altitude,  h;  the  altitude  of 
the  station,  h^;  and  the  refractivity  at  sea  level  at  the  point  in 
question,  N(0) . 

The  CRPL  Exponential  Radio  Refractivity  Atmosphere  is  presented 
(Ref.  3)  in  the  following  form: 


3-4 


PMI4C0 


WESTERN  DEVELOPMENT  LABORATORIES 


WDL-TR1993 


N(h)  =  N(h^)  exp  (h  -  (3->l) 

Pe  -  logg  [N(h^)/N(ig 

N(h_)  -  N(0)  exp  (-0.1958  h  ) 

s 

N(l)  -  N(0)  exp  (-0.1958) 


where  h 


h 

s 


N(h^) 


N(0) 


any  desired  height  (n.m.)  above  mean  sea  level, 
height  (n.m.)  of  station  above  mean  sea  level, 
refractivity  at  altitude  h,  and 

* 

sea  level  refractivity  corresponding  to  N(h  )  . 

s 


the  exponential  reference  atmosphere  is  presented  in  Fig.  3-1  for  the 
special  case  of  a  station  at  sea  level. 


The  mean  refractive  index  is  significantly  affected  by  such  vari¬ 
ables  as  season  of  the  year,  time  of  day,  and  instantaneous  weather  con¬ 
ditions  (cold  front,  warm  front,  etc.).  Although  the  theoretical  tropo¬ 
spheric  refractivity  decreases  exponentially  with  increasing  altitude, 
an  experimentally  measured  graph  or  "profile"  of  N  versus  altitude  is 
therefore  not  a  smooth  curve.  However,  the  primary  refraction  correc¬ 
tion  process  mathematically  depends  upon  the  evaluation  of  integrals 
whose  integrand  is  a  function  of  the  refractivity.  Since  integration 
is  a  smoothing  process,  the  irregularities  in  the  refractivity  profile 
tend  to  cancel  each  other;  and  the  results  are  not  very  different  from 
those  obtained  by  the  use  of  the  approximate  profile. 


Altitude  changes  in  the  refractive  index,  such  as  those  exhibited 
in  Fig.  3-l»  produce  range  and  range  rate  errors  by  creating  changes  in 
the  velocity  of  the  wave.  In  the  troposphere,  the  velocity  is  always 


*  N(hg)  can  be  estimated  from  world  charts  (Ref.  2). 
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Fig.  3«-l  CRPL  Standard  Model  for  Troposphere 

3-6 


PHILCO 


WESTERN  DEVELOPMENT  LABORATORIES 


I 

I 


I 

I 

i 


I 

I 

I 

I 


WDL-TR1993 


less  than  the  free  space  velocity  which  results  in  a  time  delay  of  the 
Signal.  In  addition,  the  increased  path  length  due  to  refraction  fur¬ 
ther  increases  the  apparent  range,  the  range  error  due  to  the  tropo¬ 
sphere  is  independent  of  frequency,  and  variations  of  the  state  of  the 
troposphere  will  cause  the  range  error  to  fluctuate  about  a  mean  value. 
The  percentage  error  due  to  tropospheric  refraction  is  the  same  for 
range  as  for  range  rate  (Doppler  frequency) . 

3.5  IONOSPHERE 

The  ionosphere  differs  from  the  troposphere  in  that  it  Contains  a 
much  higher  number  of  positive  and  negative  ions,  which  are  produced 
by  the  ultraviolet  and  particle  radiations  from  the  Sun.  The  increased 
ionization  decreases  the  effective  refractivity .  The  rotation  of  the 
earth  on  its  axis,  the  annual  course  of  the  earth  around  the  sun,  and 
the  development  of  sun  spots  all  affect  the  number  of  ions  present  in 
the  ionosphere;  and  these  in  turn  affect  the  quality  and  distance  of 
signal  transmission. 

Because  of  the  variability  of  the  ionosphere,  it  is  difficult  to 
make  exact  statements  about  the  size  of  the  tracking  error.  For  ex¬ 
ample,  at  an  altitude  of  153  n.m.,  the  ionospheric  density  has  been 
found  to  vary  as  much  as  100  percent  daily.  Throughout  a  sunspot 
cycle,  the  ion  density  has  varied  as  much  as  1600  percent  (Ref.  3). 

Different  densities  of  ionization  at  different  heights  make  the 
ionosphere  appear  to  have  layers.  The  D- layer  lies  between  heights 
of  40  and  50  n.m.  during  the  daytime.  It  is  characterized  by  low 
ionization  and  an  appreciable  absorption  effect.  The  E-layer  is  be¬ 
tween  50  and  90  n.m.  with  the  greatest  density  at  70  n.m.  and  the 
strongest  effect  in  daylight  hours.  From  90  n.m.  upward  is  the  F- 
layer  or  layers:  one  layer  at  night  and  two  layers  (Fj^  and  F2)  during 
the  daytime.  The  F-layer  has  the  highest  ionization  and,  therefore, 
the  greatest  effect  on  the  refractivity. 


3-7 


PHILCO 


WESTERN  DEVELOPMENT  LABORATOR 


WDL-TR1993 


Unlike  the  tropospheric  case,  the  effect  of  refraction  in  the 
ionosphere  is  frequency  dependent.  Both  the  angle  and  range  (range 
rate)  deviations  vary  almost  inversely  with  the  square  of  the  frequency. 
Ionospheric  refractivlty  is  a  function  of  electron  density  and  radio 
frequency.  The  ionospheric  refractive  index  in  an  Ionized  gas  (ignoring 
absorption)  follows  the  following  relationship: 


where  f^  *  signal  frequency,  and 

f^  “  plasma  frequency  or  critical  frequency. 


The  plasma  frequency,  f^,  is  a  term  frequently  used  to  express  the 
function  of  elections  and  electron  density  in  determining  the  index  of 
refraction.  The  actual  relationship  is  given  as 


f 

c 


(3-3) 


where  =  electron  density, 

e  s  electron  charge,  and 
m  o  electron  mass. 


Substituting  values  for  e,  m,  and  n  in  equation  (3-3)  and  combining 
with  equation  (3-2),  the  following  formula  is  obtained: 


n 


(3-4) 


The  electron  density,  N^,  must  be  known  in  order  to  calculate  the 
refractive  index  by  equation  (3-4).  Electron  density  profiles,  giving 
for  various  altitudes  and  geographic  locations,  can  be  obtained  by 
ionosphere  soundings;  that  is,  actual  measurements  via  rockets  or  probes. 
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In  the  absence  of  measured  electron  density  profiles,  the  Chapman  for¬ 
mula  may  be  used  as  a  model  for  the  electron  density  in  the  ionosphere. 
The  Chapman  formula  is  the  following: 


N  ■  N  exp  -i  (1  -  Z  -  e^  see  5) 
e  m  2  ’ 


m 


H 


H  = 


kl 

mg 


where  N 

m 

h 

h 

m 

k 

T 

m 

g 


maximum  electron  density  when  the  sun  is  overhead, 
solar  zenith  angle, 

any  desired  height  in  nautical  miles  above  mean  sea  level, 

height  of  maximum  density  N  , 

in 

Boltzman's  constant, 
absolute  temperature  (at  altitude  h), 
molecular  mass  of  the  gas,  and 
gravitational  constant. 


The  average  variations  of  the  critical  frequencies  of  the  iono¬ 
spheric  layers  are  sufficiently  well  known  to  permit  long-range  pre¬ 
dictions  (accurate  to  within  about  25  percent)  of  average  conditions 
expected  on  "quiet"  days.  Such  predictions  are  published  three  months 
in  advance  by  the  National  Bureau  of  Standards.  However,  corresponding 

values  of  h  and  H  must  be  known  to  obtain  an  electron  density  profile 
m 

even  when  using  the  Chapman  formula.  Table  3-1  indicates  the  values  of 
the  parameters  in  Chapman's  formula  obtained  from  extensive  pulse 
soundings  and  statistical  analyses  during  periods  between  1948  and  1954 
at  Washington  D.  C>  (Ref>  4,5).  Figure  3-2  shows  some  ionospheric 
electron-density  profiles  computed  from  Chapman's  theory.  The  three 
solid  curves  represent  the  2.5,  50.0,  and  97.5  percent  confidence  limits 
of  the  statistical  analysis.  The  broken  curves  are  the  noon  and  mid¬ 
night  results  of  a  single  experiment  measuring  refractivity.  The 
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critical  (plasma)  frequency  corresponding  to  any  point  on  a  curve  may  be 
read  on  the  upper  horizontal  scale,  while  the  lower  horizontal  scale 
represents  the  electron  density  corresponding  to  each  critical  frequency. 
Figure  3-3  was  obtained  by  using  the  information  from  Fig.  3-2  in  equa¬ 
tion  (3-4)  to  obtain  a  refractivity  profile  for  an  operating  frequency 
of  400  Me. 


TABLE  3-1 


MAJOR  CHARACTERISTICS  OF  IONOSPHERIC  LAYERS 


Layer 

h 

m 

average 
height 
(n.  miles) 

H 

scale 
height 
(n.  miles) 

Critical  frequency  and 

Maximum  electron  density” 

2.5 

percentile 

rms  or 
average 

97.5 

percentile 

0.7 

1.1 

1.3 

D 

46 

2.4 

6.0 

15 

20 

0.9 

3.3 

4.8 

E 

58 

5.9 

9.2 

134 

290 

2.7 

4.4 

5.0 

^1 

100 

16 

90 

240 

310 

2.7 

5.3 

10.6 

^2 

150 

38 

90 

354 

1400 

*  Upper  value  is  critical  frequency  in  megacycles  per  second;  lower 
value  is  maximum  electron  density  in  units  of  10^  per  cubic  meter. 
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3 . 6  STANDARD  CORRECTIONS 

Corrections  for  tropospheric  angle  refraction  can  be  made  in  sev¬ 
eral  ways.  One  method,  based  on  the  surface  index  of  refraction,  Claims 
prediction  accuracies  of  0.02  milliradian  at  a  5-degree  antenna  eleva¬ 
tion  angle.  A  simpler  model,  using  the  yearly  average  of  the  station 
Surface  index  of  refraction,  also  reduced  the  refraction  error  signifi¬ 
cantly.  Correction  techniques  using  the  world-wide,  all-year  average 
are  also  capable  of  reducing  the  error,  although  not  as  much  as  the 
above  methods.  The  most  practical  correction  method  for  uSe  in  tracking 
vehicles  seems  to  be  the  yearly  average  of  the  station  surface  index  of 
refraction. 


Compensation  for  ionospheric  effects  can  be  carried  out  by  methods 
similar  to  those  used  for  tropospheric  effects.  By  employing  frequent 
ionospheric  soundings  to  determine  the  profile  of  the  ionosphere,  the 
amount  of  expected  refraction  can  be  calculated.  A  less  accurate  method 
would  be  to  use  typical  profiles  either  from  past  soundings  at  the  sta¬ 
tion  or  from  other  stations  near  the  same  latitude.  If  this  type  of 
correction  is  used,  a  day  profile  and  a  night  profile  are  needed.  Using 
frequent  soundings,  the  range  rate  and  angle  error  can  be  reduced  by 
as  much  as  80  percent.  Since  the  ionization  of  the  ionosphere  can  vary 
over  such  wide  ranges,  a  correction  made  from  average  profiles  does  not 
seem  to  be  very  reliable. 

Both  the  CRPL  model  for  the  troposphere  and  the  Chapman  model  for 
the  ionosphere  are  attempts  to  establish  a  "standard"  atmosphere,  al¬ 
though  no  such  thing  exits.  Not  only  are  there  changes  from  day  to  day, 
but  there  are  also  turbulence  layers  as  well  as  occasional  clouds  and 
precipitation  particles.  The  use  of  a  standard  atmosphere  model  does 
not  give  perfect  results;  however,  continuing  investigations  are 
yielding  improved  solutions  to  the  problem. 


3-13 


PHILCO 


WESTERN  DEVELOPMENT  LABORATORIES 


WDL^TR1993 


SECTION  4 
REFRACTION 


4.1  GENERAL 

The  purpose  for  developing  such  models  aS  the  CRPL  Exponential 
Atmosphere  and  the  Chapman  formula  for  electron  density  in  the  iono¬ 
sphere  is  to  enable  a  reconstruction  of  the  signal  path  by  the  laws 
of  refraction.  Refraction  correction  formulas  for  optical  or  radio 
tracking  instrumentation  are  generally  based  on  the  use  of  geometrical 
optics  in  a  spherically-stratified  terrestrial  atmosphere. 

4.2  GEOMETRICAL  OPTICS 

Geometrical  optics  is  the  limiting  case  of  physical  optics  (based 
on  Maxwell's  electromagnetic  wave  equation)  as  the  wavelength  approaches 
zero,  and  its  use  is  valid  in  ray  tracking  problems  where  the  refractive 
index  does  not  vary  much  in  a  distance  comparable  to  the  wavelength. 

The  radio  wavelength  is  about  10  feet  at  100  Me  and  about  1  foot  at 
1,000  Me,  while  the  atmospheric  index  at  these  frequencies  normally 
changes  only  by  amounts  on  the  order  of  parts,  or  tenths  of  parts,  per 
million  over  such  distances.  Hence,  the  use  of  geometrical  optics  in 
tracking  problems  appears  well  founded. 

At  the  interface  of  two  optical  media  with  differing  indexes  of 
refraction,  n^^  and  n2^  the  direction  of  a  straight-line  ray  path  will 
change.  The  incident  ray,  the  line  normal  to  the  interface,  and  the 
refracted  ray  are  all  coplanar;  if  9^  and  angles  between 

the  normal  and  the  two  rays,  then  Snell's  Law  states  that 

sin  ■  n2  sin 

If  a  refractive  medium  has  an  index  of  refraction,  n,  which  de¬ 
pends  only  on  the  distance,  r,  from  a  given  point  (spherical  sytmnetry), 
then  Bpuguer's  Law  states  that  the  ray  path  is  described  by 
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n  slti  ^  constant 
r 

where  n  and  r  pertain  to  any  point  on  the  ray  and  ^  denotes  the  angle 
between  the  radius  vector  and  the  tangent  to  the  ray  path  at  the  given 
point  (((  is  the  GOmplement  of  angle  in  Fig.  4«1). 


4.3  REFRACTION  GEOMETRY 

Figure  4-1  shows  the  typical  spherical-earth,  ray  refraction  geo¬ 
metry  in  the  plane  determined  by  the  center  of  the  earth  0,  a  ground 
station  S,  and  the  ray  path  ST. 


Let 


R 

h 


Q1TQ2 

SQ3 

STQ, 


E 


E^+AE 

T 

T 

7 


T 


Ni(hs) 

N<hj.) 

N(h^) 


s  mean  radius  of  the  earth, 

>  height  of  station  above  mean  sea  level, 
height  of  target  above  mean  sea  level, 

=  line  tangent  to  the  ray  at  T, 

B  line  tangent  to  the  ray  at  S, 

■  straight  line  through  S  and  T  (radial  line), 

>  observed  elevation  angle, 

«  true  elevation  angle, 

*  elevation  angle  error,  »  E^  -  E  >  0, 

declination  angle  as  observed  from  satellite, 
c  declination  angle  as  observed  from  station, 

•=  declination  angle  error, 

=  total  bending  angle,  7  » 

=  refractivlty  at  the  ground  station, 
s  refractivlty  at  the  target,  and 

s  refractivlty  at  an  altitude  h.  above  mean  sea  level. 


If  the  refractivlty  profile  for  the  function  N(h)  is  not  known  exactly, 
it  may  be  approximated  by  a  standard  atmosphere  as  in  Section  3. 
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4.4  TOTAL  BENDING  AND  ELEVATION  ANGLE  CORRECTION 

A  method  is  now  presented  for  computing  the  bending  of  radio  rays, 
7,  through  an  atmosphere  of  "known”  refractive  index  distribution,  under 
the  usual  assumption  that  surfaces  of  equal  index  are  spherical  and  con¬ 
centric  with  the  earth. 


Consider  a  ray  entering  an  infinitesimal  layer  of  thickness  dh  at 
an  angle  ^  (Fig.  4-2).  If  K  is  the  radius  of  curvature  of  the  refracting 
ray^  then  (Ref.  12): 


1 

R 


1  dn 

-  ■—  cos  3 
n  dh 


(4-1) 


where  n  is  the  index  of  refraction.  The  length  of  the  ray  path  in  the 
layer  is 


K  dy  «  CSC  3  dh. 


(4-2) 


This  may  be  seen  by  considering  the  arc  of  a  circle  and  two  tangents 
intersecting  each  other  in  an  angle  dy.  The  angle  formed  at  the  center 
of  curvature  by  the  two  radii  must  also  be  dy  since  all  angles  must 
total  2n.  The  length  of  the  enclosed  arc  is  then  K  dy . 

If  the  layer  is  thin,  the  arcs  may  be  approximated  by  the  tangent 
lines;  and  the  triangle  with  a  hypotenuse  of  length  dh  esc  3  1^  formed. 
The  hypotenuse  length  is  very  nearly  the  same  as  the  length  of  the  en¬ 
closed  arc.  Therefore,  equation  (4-2)  is  established. 


dy  ■  —  P 

n  dh 


(4-3) 


The  dy's  of  all  elementary  layers  are  directly  additive, 
in  Fig.  4-2.  It  follows,  therefore,  that  the  contribution  to 
bending,  y,  due  to  a  layer  bounded  by  the  heights  end 


as  shown 
the  total 
is 
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’'k,k+l 


A  ^ 

n  dh 


cot  dh  . 


(4-4) 


A  convenient  representation  for  the  angle  Pj^  is  now  required.  If 
the  ray  departs  from  the  earth's  surface  at  the  elevation  angel  E^, 
Bouguer's  Law  for  spherical  stratification  states: 

n^(R  +  h^)  sin  (90  «  E^)  *  ’^(^  (^0  “  “  constant 

"sCR  +  h^)  cos  E^  •  nj^(R  +  hj^)  Cos  pj^ 


cos  E 


cos  £ 


CO.  ^ 


n.  TR  +  hn  r  -I  'r  + 

T  TTT  [1  -  NCh^)  -  N(h,)  10-®  — 

8  L  sj  *-  L 


R  +  h 

s 


(4-5) 


Having  the  angle  Pj^,  the  integral  equation  (4-4)  can  be  evaluated.. 
Weisbrod  and  Anderson  (Ref.  12)  have  shown  that  the  integral  has  the 
following  value  for  each  interval  h^^  to 

^  l»(Wl-  l”<Vl 

'k,k+l  500  (tan  +  Can  Pj^) 


When  these  small  bending  angles  are  summed,  the  total  bending  angle  is 
obtained: 


_  l»(Wl  - 

^  ^  500  (tan  Pj^^^  +  tan  p^) 


(4-7) 
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The  total  ray  bending  in  a  spherically--stratified  atmosphere  can 
be  evaluated  by  approximating  the  refractivity  profile  by  a  series  of 
straight-line  segments,  computing  the  bending  due  to  each  of  the  cbr- 
responding  Spherical  shells,  and  adding  the  results  to  obtain  the  total 
bending. 


The  first  step  is  to  select  a  series  of  significant  altitudes 

h»,  h,  ,  h.,,  h.,,  ...  h .  (h  =  h  ,  altitude  of  station:  h.  =  h^, 

0  12  3  JOS  j  t 

altitude  of  target)  such  that  between  each  consecutive  pair  the  re- 
fractivity  profile  is  well  approximated  by  a  straight  line.  Usually, 

ie 

a  number  of  height  levels  between  50  and  60  will  suffice  for  the  total 
atmospheric  height  encountered  in  missile  and  satellite  work.  This  will 
yeild  sufficient  accuracy  while  retaining  the  characteristics  of  the 
relatively  finer  structure  of  the  atmosphere. 

The  next  step  is  to  tabulate  the  refractivity  and  calculate  the 
angles  as  in  equation  (4-5),  then  to  use  this  information  to  cal¬ 
culate  the  total  bending  angle  7.  It  is  known  that: 

7  = 

Having  already  obtained  a  solution  for  7,  solution  for  will  give  an 
answer  for  AE.  Referring  again  to  Fig.  4-1,  the  law  of  a  sines  gives; 

R  +  h  R  +  h 

_ t  _  _ s 

cos  (E^  -  AE)  ^  cos  (E^  +  AE^) 


*  The  effect  of  varying  the  height  increment  used  for  refraction^ 
correction  calculations  is  discussed  in  Section  6. 
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R  +  h 


+  AE^  »  arc  cos  cos  (E^  -  AE) 


(4-8) 


By  Bouguer's  Law, 


(R  +  h^)  cos  E^  *  (R  +  h.)  n  cos  E  , 
t'  t  T  s'  s  o’ 


E^  =  are  cos 


n  R  + 

^8  _ 

R  + 


Cos  E 
'■t 


Hence,  combining  equations  (4-8)  and  (4-9): 


R  +  h  n  R  +  h 

AE_  “  arc  cos  r— -  cos  (E  -  AE)  -  arc  cos  — ^  ^  cos  E  .  (4*10) 

i'  K+il  O  nK+il  o 

t  t  t 

Expanding  the  trigonometric  functions  of  this  expression  by  series  and 
including  seventh-power  terms  in  the  cosines >  as  Is  necessary  to  retain 
required  accuracy: 

^  Fr  +  h  n  ~[  \  2j-l  /  cos  (E  -  AE)\  2j-l 


(4*9) 


R  +  h  n 

im;  s;  ® 


(1  +  AE  tan  E  ) 
n  y  o 


■4 


(4-11) 


PHIICO 


WESTERN  DEVEUPPMENT  LABORATORIES 


WDL-TR1993 


Now  since 


+  [N(hp  -  N(hg)jlG"j  ^  I  +  (2j-l)jN(hj.)  -  N(h_)JlO' 


(4-12) 


and 


(I  +  AE  tan 


2j-l 


1  +  (2j-l)  AE  tan  E^, 


(4-13) 


Gombining  equations  (4-12)  and  (4-13)  results  in 


'nt\2  j-1 

(1  +  AE  tan  E 
n  y  '  o 

S'' 


)2j"^  ^  1  +  (2j-l)|jN(hj.)  -  N(h^g  10'^ 


+  AE  tan  I  (4-14) 


•1 


Equation  (4-12)  for  AF.^  thus  becomes: 


-  ^  ^  tan  E 


(2j-l) 


rR  +  h  n 

g-'  V  ^  ~  cos  E 
R  +  n^  o 


t  t 


(4-15) 


Since  AE  =  7  -  AE^,  substitution  of  the  last  equation  yields: 


7  + 


|^(hj.)  -  N(h^^  10 


AE 


1) 


R  +  h  n 
s  s 

- ^  cos  E 

R  +  h  n  o 


-l2j.l 


1  -  tan  E 


4  r- 

ZIR  +  h  n 

j=l  U  t  C 


cos  E 


2j-l 


(4-16) 
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Furthermore: 


I 


Cj(2j-1) 


n  R  +  h 
8  8  ^ 
—  r’-'V '  'r  -  cos  E 
n  R  +  h  o 


^2j-l 


J-1 


(2j-l) 


1  + 


(2j-l)  [N(h^)  -  N(h^3  10 


-6 


R  +  h 


s 


a  .  ,  cos  E 
R  +  o 


2j-l 


thus  equation  (4-16)  becomes: 


^  7  *  [M(h^)  -  N(h^)1 10-^ 


1  -  tan  E  Y- 
o  / 


(4-17) 


where 


^1  “  Z  Cj(2j-1)  X 

J-1 


2j-l 


Y2  ^  +  [N(h^)  -  N(hj.)]l0’^  ^  Cj(2j-1)  X 

j-1 


2  2j-l 


X  - 


R  +  h 

s 

R  + 


cos  E 

o 


^  1,  1/6,  C3  -  3/40,  5/112 

By  using  the  above  equations,  it  is  possible  to  reconstruct  any 
ray  path  and  to  obtain  a  correction,  for  the  refracted  elevation 
angle . 
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SECTION  5 

DOPPLER  CORRECTIONS 

5.1  ELEVATION  ANGLE  REFRACTION  ERROR 

Let  T  represent  a  satellite  transmitter  location  and  the  magni¬ 
tude  of  its  velocity  projected  onto  the  plane  determined  by  the  ray 
path.  The  magnitudes  of  the  components  of  the  satellite  velocity  in 
the  various  directions  shown  in  Fig.  5-1  are 


measured  radial  velocity 

.  «.  V 

r 

=  V 

P 

cos  (T|f  +  AE^) 

(5-1) 

true  radial  velocity 

=  V 

P 

00  S  i|r 

(5-2) 

apparent  radial  velocity 

=  V 

a 

=  V 

p 

cos  (ijr  -  AE) 

(5-3) 

The  quantity  of  interest  is  the  error  in  the  true  velocity  component. 
The  error  is  the  difference  between  the  true  velocity  component  and  the 
measured  velocity  component: 

-  V^.  (5-4) 

Substituting  values  from  equations  (5-1),  and  (5-2)  into  equation  (5-4) 
gives: 


AVj  =  Vp  cos  t  ~  Vp  cos  (i|f  +  AE.J,) 

=  V  cos  -  V  cos  *  cos  AE_  +  V  sin  W  sin  AE_ 
p  p  T  p  T 

=  Vp  cos  ^  (1  -  cos  M.p)  +  Vp  sin  i|r  sin  AE^.  (5-5) 


Since  AE^  is  a  very  small  angle,  the  trigonometric  functions  can  be 
expanded  by  series  to  the  second  order: 


2 

AVd  ^  Vp  0E^  sin  \|f  +  •|(AE.j.)  cos 


(5-6) 
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STATION  COORDINATE  SYSTEM  CONFIGURATION  OF  VELOCITY 

VECTORS  AT  THE  SATELLITE 


/ 

/ 

/ 


FiR,  5»I  Geometry  for  Doppler  Refraction  Error  Analysis 
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For  cases  where  >|r  is  reasonably  different  from  zero,  this  expression 
may  be  reduced  to 


AV^  ^  Vp  AE^  sin  (5-7) 

5.2  ONE-WAY  DOPPLER  REFRACTION  ERROR 

From  Section  2,  the  one-way  Doppler  frequency  shift  is 

f  . 

£  o  ds 

d  c  '^t  dt  ‘ 

The  time  rate  of  change  of  the  signal  path  length  is  ds/dt.  This  signal 
path  length  is  not  too  different  from  the  slant  range,  R;  therefore,  the 
time  rate  of  change  of  both  s  and  R  should  not  be  appreciably  different 
and  could  be  approximated  as 

ii  is  =  R 

dt  dt 

By  also  substituting  R  =  V  cos  0,  where  0  is  the  angle  between  the  true 
velocity,  V,  and  the  position  vector  of  the  satellite  relative  to  the 
station,  the  resulting  formula  for  one-way  Doppler  frequency  shift 
becomes 


f  =  (f  /c)  n  V  cos  0 
dot 

(f^/c)  n^. 


(f  /c)  n^  V  cos  * 
o  t  p  ’ 

*  (f^/c)  n^.  V  cos  0  cos  i|f 


(5-8) 
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where  0 


f 


d 


£ 

o 

c 


angle  between  V  and  the  ray  path  plane, 

Doppler  shift,  or  difference  between  the  transmitted  fre¬ 
quency  and  the  apparent  received  frequency, 

transmitted  frequency, 

velocity  of  light  in  a  vacuum,  and 

index  of  refraction  at  the  target  location  corresponding 
to  the  transmitted  frequency  f^. 


Differentiating  equation  (5-8),  the  error  in  Doppler  frequency  due  to 
the  velocity  error  is 


Af ,  =  (f  /c)  n  cos  f  AV  (5-9) 

dot  P 

where  f^,  c,  n^,  and  are  assumed  to  be  errorless  quantities. 

By  differentiating  equation  (5-2),  the  resulting  expression  is 


AV .  »  AV  cos  if 

dp 


Equation  (3-9)  becomes 


^d  "  ^  "t  ^'^d 


which,  when  combined  with  equation  (5-6),  yields: 

"t  ~  ”t  ^  10. 

Calculation  of  the  angle  has  been  discussed  in  Section  4.  With 
this  result,  it  is  possible  to  calculate  the  Doppler  frequency  error 
due  to  refraction  in  the  case  of  one-way  Doppler  transmission. 
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5.3  TWO-WAY  DOPPLER  REFRACTION  ERROR 

Using  the  two-way  Doppler  expression,  equation  (2-6),  the  Doppler 
frequeney  error  due  to  refraction  in  two-way  transmission  is 


5.4  DOPPLER  ERROR  CORRECTION  EXPRESSIONS 

The  Computation  of  errors  induced  by  refraction  of  dOppler  trans 
missions,  as  developed  in  previous  sections,  requires  a  method  for  cal¬ 
culating  the  angle  \lr.  With  reference  to  Fig.  5-1,  construct  the  Station 
centered,  left-handed  Cartesian  coordinate  system  (Sj^,  S^,  S^)  and  Cor¬ 
responding  unit  vectors  (Sj^,  §2,  S^)  respectively,  where  the  S^  axis  is 
directed  toward  true  north,  the  $2  axis  east,  and  the  S^  axis  vertical 
so  that  the  vector  S^  points  along  the  vertical  axis  and  is  given  by 
the  vector  product  =  §2  «>S^. 

For  any  given  satellite  pass  at  time  tj^,  the  following  observed 
parameters  are  defined  as: 

Rj^  =  slant  range  (range  from  station  to  vehicle), 

Aj^  *  azimuth  angle  (measured  from  true  north),  and 

Ej^  =  elevation  angle  (measured  from  the  (Sj^,  S2)  plane). 

The  position  vector,  P^,  from  the  station  to  the  vehicle  at  time  is 
expressed  as: 

Pj^  =  (Rj^  cos  cos  Ej^)  Sj^  +  (Rj^  sin  cos  Ej^)  §2  +  (Rj^  sin  Ej^)  S^-  (5-12) 

Furthermore,  the  azimuth  vector,  A^^,  the  normalized  projection  of 
vector  Pj^  on  the  (S^,  S2)  plane,  is  defined  as: 


Aj^  -  (cos  A^)  +  (sin  A^)  §2. 
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the  curved  ray  path,  S,  of  the  doppler  transmission  is  taken  to  be  in 
the  plane  formed  by  the  vectors  and  so  that  the  vector  a 
unit  normal  vector  to  this  plane,  is  given  as: 


N, 


■nrr  ■ 


Sin  E. 


k  (5-14) 


or 


-  (Sin  A^)  -  (cos  Aj^)  S^. 


(5-lS) 


The  vehicle  velocity  vector,  Vj^,  is  expressed  by  a  finite  vector 
difference  between  two  consecutive  position  vectors  as:: 


(P, 


k+1 


^k> 


(5-16) 


The  component  of  the  velocity  vector,  Vj^,  in  the  ray  path  plane  is  rep¬ 
resented  by  the  vector  V^,  which  is  formed  by  the  vector  product  of 

V,  and  N.  so  that: 
k  k 


V.  (g) 


(5-17) 


The  angle  \|f  is  defined  as  the  angle  between  the  vector  and  the  vector 

V, .  The  vector  V,  is  defined  as  the  true  (errorless)  statlon-to-vehicle 
Q  O 

velocity  vector, 
affect  only  elevation  angles. 


lies  in  the  ray  path  plane  since  refraction  errors 


The  angle  t  expressed  as: 


y  ■  arc  cos 


rv;; 


[^df 


(5-18) 
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However,  in  the  normal  case  of  data  reduction,  data  is  not  available 
for  the  vector  V^.  The  vector  can  be  approximated  (in  direction,  only) 
by  normalizing  the  vector  corrected  for  elevation  errors  due  to  re¬ 
fraction.  (Approximate  values  for  the  elevation  errors,  are  found 
in  Fig.  5-6.) 

Denote  the  corrected  vector  Pj^,  which  is  normalized  to  magnitude 
one,  as  the  vector  so  that, 

Sj  +  cos  (E^  +  §2 

+  Sin  (Ej^  +  AE^)  §3.  (5-19) 

Then  the  angle  \lf  is  expressed  as: 


i),  «  C( 


cos  A,  cos  (E  +  AE 


^  = 


(5-20) 


5.5  ORBITAL  VELOCITY  ERROR 

The  magnitude  of  the  error  induced  in  orbital  velocity  due  to  re¬ 
fraction  effects  is  denoted  by  AV. 


The  values  of  the  angle  \|r,  obtained  from  equations  (5-10)  and  (5-11), 
determine  the  results  for  Af^,  so  that  AV^  is  given  by  AV^  =  \ 
where  \  is  the  wave  length  corresponding  to  f^. 

From  equation  (5-2)  and  the  geometry  in  Fig.  5-1,  the  errors  in 
magnitude  of  vectors  and  V  are  given  as 

AV  AV ,  sec  \lf 

”  p  -  d 

AV  »  AV  sec  i 

p 
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where  i  is  the  angle 


between  V  and  V 

P 


and  is  given  by: 


cos 


(5-23) 


Combining  the  above  results,  the  error  in  orbital  Velocity  becomes 


AV  * 


XAf . . 

d 


(5-24) 


Figure  5-2  and  Table  5-1  give  examples  of  results  of  accuracy 
expected  for  Doppler  errors  based  on  the  development  of  the  above  and 
previous  analyses. 


5.6  REFRACTION  CORRECTION  FOR  ORBITAL  ERROR 

In  determining  the  correction  for  refraction  effects,  let 

f^^(t)  ■  observed  one-way  Doppler  frequency  shift  at  time  t, 

Ep(t)  =  observed  elevation  angle  at  time  t, 

f^(t)  =  Doppler  frequency  shift  at  time  t  in  absence  of 

atmosphere, 

E(t)  B  elevation  angle  at  time  t  in  absence  of  atmosphere, 
and 

M  »  refraction  error  in  apparent  elevation  angle  E 
at  time  t.  ° 

It  is  obvious  that: 

E^(t)  »  E(t)  +  m. 


If  the  refraction  error  is 


AE 


E'(t)At, 
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then 


E^(t)  -  E<t)  +  E'(t)At  «  E(t  +  At) 


or,  equivalently: 


E(t)  *  E^(t  -  At)  -  E^(t)  ^  E^(t)At, 


where 


At  » 


IaeI 

E '(t)  * 


(5‘25) 


Equation  (5-25)  shows  how  to  recalibrate  the  time  scale  so  that 

the  true  (vacuum)  value  of  the  elevation  angle,  E,  may  be  expressed  as 

the  value  of  the  observed  elevation,  E  ,  at  some  other  time.  It  follows 

’  o* 

that  the  same  transformation  of  the  time  scale  will  give  the  true  value 
of  the  Doppler  shift,  f^,  in  terms  of  the  observed  Doppler  shift,  f^^, 
at  the  transformed  time  value: 

fj(t)  •=  f^^(t  -  At).  (5-26) 


Figure  5-3  shows  the  application  of  this  method.  The  observed 
values  f^^  and  E^  are  represented  by  solid  curves, while  f^  and  E  are 
represented  by  dashed  lines. 

The  corrected  elevation  curve  was  obtained  from  the  refraction 
formulas,  equations  (4-5)  and  (4-7).  The  corrected  Doppler  shift  curve 
was  plotted  by  taking  time  shifts  (At)  corresponding  to  the  time  dif¬ 
ference  between  the  two  elevation  curves.  This  same  time  difference, 
when  applied  to  the  observed  Doppler  curve,  forms  the  corrected  fre¬ 
quency  shift  curve  shown  in  Fig.  5-3. 
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SECTION  6 

EXPECTED  ACCURACY  OF  DOPPLER  MEASUREMENTS 


6 . 1  GENERAL 

Determination  of  accuracies  in  any  refraction  correction  procedure 
is  considered  one  of  the  important  aspects  Of  the  present  study.  Com¬ 
putations  were  carried  out  (Ref.  3)  to  demonstrate  the  effects  of  the 
size  of  increments  in  the  expression  used  to  approximate  typical  refrac¬ 
tion-correction  integrals.  Similar  calculations  were  made  to  illustrate 
the  effect  of  rather  drastic  approximations  to  the  shape  of  ionospheric 
layers  as  well  as  to  determine  the  effect  of  changing  the  height  of 
Such  layers.  Based  upon  these  results,  an  analysis  using  results  of 
numerical  work  by  Counter  and  Reidel  (Ref.  5)  was  conducted  to  estimate 
accuracies  expected  of  refraction-correction  procedures.  A  determina¬ 
tion  of  the  errors  resulting  from  uncertainty  in  the  local  index  at  the 
target  was  the  basis  of  another  analysis. 

6.2  EFFECT  OF  INCREMENTAL  ALTITUDE  CALCULATIONS 

In  obtaining  the  correction,  AE,  to  the  apparent  elevation  angle, 
E^,  the  total  bending,  7,  must  be  calculated  through  use  of  a  finite- 
summation  approximation  to  an  integral.  Often  the  values  of  h^  in  this 
expression  are  conveniently  chosen  to  be  equally  spaced  with  a  certain 
altitude  increment.  For  example,  Weisbrod  and  Anderson  (Ref.  12)  state 
that  the  resultant  error  in  tropospheric  bending  is  only  about  9.6  per¬ 
cent  for  a  lO-mile  increment. 

An  analysis  was  conducted  (Ref.  3)  to  check  the  effect  of  selected 
increments  upon  ionospheric  refraction  corrections.  The  results  of 
these  calculations  indicate  that  the  increment  chosen  for  calculation 
has  a  smaller  effect  upon  the  results  than  might  be  expected  at  first. 
However,  wherever  reasonably  high  accuracies  are  required  in  the  cor¬ 
rection  procedures,  smaller  increments  should  be  used  to  preserve  the 
details  of  the  refractivity  profile,  assuming  the  profile  itself  is 
sufficiently  well  known  to  warrant  this  procedure* 
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6.3  CHAPMAN  MODEL  LINEAR  APPROXIMATION  ERRORS 

The  effect  of  ionospheric  layer  shape  was  investigated  (Ref.  3) 
by  approximating  the  refractivity  profile  corresponding  to  the  97.5 
percentile  curve  of  Fig.  3^3  by  a  relatively  small  number  Of  straight 
line  segments  (in  this  case,  five).  The  approximate  profile  thus  ob¬ 
tained  was  compared  to  the  original  profile  in  Fig.  6-1.  the  linear 
approximation  to  the  refractivity  profile  is  given  by  the  following 
function  (for  400  Me) : 


Altitude 

(n-m* ) 

Approximate  Refractivity  (N  units) 

150  <  h 

<  600 

-1550  exp  (-.0096  h) 

97  <  h 

<  150 

-71  exp  (+.011  h) 

76  <  h 

<  97 

-.027  exp  (+0.092  h) 

^  “  N(h) 

53  <  h 

<  76 

-601  exp  (-.039  h) 

45  <  h 

<  53 

-45  X  lO'®  exp  (+0.37  h) 

In  a  previous  experiment  (Ref.  3),  it  was  observed  that  AE  cal¬ 
culated  from  the  true  profile  differed  from  AE  calculated  from  the 
linear  approximation  by  less  than  0.05  mr  in  most  cases.  The  linear 
approximation  is  sufficiently  accurate,  therefore,  to  use  in  the  cal¬ 
culation  of  bending  and  elevation  angle  error. 

6.4  EFFECT  OF  IONOSPHERIC  LAYER  HEIGHT  CHANGES 

The  effect  of  changes  in  the  height  of  the  maximum  chosen  for  the 
F-layer  of  the  ionosphere  was  investigated  (Ref.  3)  for  f^  ■  100  Me 
by  using  the  linear  approximation  profile  for  comparison  purposes. 

The  F-layer  portion  of  this  profile  was  moved  up  in  height  by  10  and 
25  n.m.,  and  the  resulting  elevation  angle  errors  were  calculated. 


It  was  concluded  that  the  choice  of  the  height  of  the  ionspheric 
layer  has  a  most  important  effect  upon  the  values  attained  for  refrac¬ 
tion  corrections. 
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6.5  EFFECT  OF  ANGULAR  ERROR  STANDARD  CORRECTIONS 

In  addition  to  the  foregoing  analyses,  expected  accuracies  of 
standard  corrections  were  investigated  with  the  aid  of  numerical  work 
by  Counter  and  Riedel  (Ref.  5).  For  example,  Figs.  6-2  through  6-5 
summarize  the  results  of  a  standard  refraction  correction  applied  in 
a  case  where  f^  =  400  Me.  The  graphs  Show  the  expected  value  of  the 
elevation  angle  deviation,  E^  -  E,  the  standard  correction,  AE,  applied, 
and  the  standard  deviation  of  the  bias  remaining  in  the  elevation  angle, 
with  and  without  the  application  of  a  correction  based  on  a  measured 
value  of  the  index  of  refraction  at  the  ground  at  the  time  of  the  ob¬ 
servation  of  the  target.  All  Of  the  curves  are  drawn  as  contours  of 
constant,  true  elevation  angle,  E. 

The  fact  that  the  curves  in  Figs.  6-4  and  6-5  are  not  all  zero  for 
an  arbitrary  altitude  indicates  the  failure  Of  the  standard  atmosphere 
to  conform  precisely  to  the  actual  atmosphere.  From  Fig.  6-4,  it  may 
be  seen  that  the  angular  refraction  corrections  for  the  tropospheric 
(h  <  40  n.m. )  bias  errors  should  have  good  accuracy  (better  than  0.05 
mils)  at  elevation  angles  greater  than  30  degrees.  It  may  be  concluded 
from  Fig.  6-5  that  the  added  knowledge  of  surface  index  of  refraction 
should  give  good  accuracies  for  elevation  angles  greater  than  6  degrees. 

For  ionospheric  effects  on  angle  measurement,  however,  the  standard 
correction  will  not  give  comparable  results.  Data  from  ionospheric 
soundings  will  be  necessary  to  obtain  the  greater  accuracies  which  are 
necessary  for  precise  Doppler  velocity  measurements. 


6.6  EFFECT  OF  DOPPLER  MEASUREMENT  STANDARD  CORRECTIONS 

Referring  to  equation  (5-10),  it  can  be  seen  that  the  maximum 
must  be  on  the  order  of  10  mils  or  less  if  the  accuracy  in  the  velocity 
correction  is  to  approach  0.01  foot  per  second.  However,  Figs.  6-4  and 
6-5  show  that  it  is  doubtful  whether  this  accuracy  can  be  achieved  under 
the  best  conditions.  Even  if  the  accuracy  is  relaxed  to  0.1  foot  per 
second,  a  standard  correction  for  tropospheric  effects  will  not  yield 
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the  desired  accuracy  without  measuring  the  local  surface  index  or  re* 
fraction.  An  accuracy  of  0.1  foot  per  second  is  suitable  (according 
to  Ref.  3)  for  ionospheric  predictions  and  other  similar  data,  but  only 
for  angles  greater  than  about  15  degrees. 


6.7  EFFECT  OF  LOCAL  INDEX  UNCERTAINTY  AT  TARGET 

An  additional  error  is  introduced  if  the  effect  of  the  local  index 
at  the  target  is  not  considered  when  converting  from  Doppler  frequency 
to  Doppler  velocity.  From  equation  (5-10),  it  is  evident  that  the 
effect  of  the  local  index  can  be  important  in  making  the  conversion. 
Assuming  a  target  in  the  ionosphere  transmitting  to  a  ground  station, 
the  ionospheric  index  from  equation  (3-4)  is 


n^.  »  -Jl  -  80.5  (Ng/f^)  isi  1  -  40.2  (N^fJ)  (6-1) 


or,  from  equation  (3-2) : 


(6-2) 


Now  since  R  =  f A,  where  \  -  c/f  n^, 
d  ’  /  o  t’ 


R 


^o"t' 


(6-3) 


Differentiating  this  expression  to  find  the  dispersion  corresponding 
to  an  uncertainty  in  knowledge  of  local  index: 


fjC 

^  =  7-2  ^t* 
Vt 
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Therefore: 


M 

A 


(6-4) 


in  the  ionosphere,  where  this  effect  can  become  important^  the  expres¬ 
sions  for  ionospheric  index,  (6-1)  and  (6-2),  yield  by  differentiation: 


40.5 


M 

e 


Combining  equations  (6-4)  and  (6-5): 


(6-5) 


•  •  c  •  40  5 

AR«-R-|Af  «-R  M  .  (6-6) 

f.i  c  r  ® 

o  o 

Hence,  the  velocity  error  will  depend  upon  the  uncertainty  in  knowledge 
of  the  critical  frequency  or  electron  density.  Figure  6-6  shows  these 
velocity  errors  as  a  function  of  operating  frequency  for  different 
values  of  uncertainty  in  critical  frequency.  The  maximum,  mean,  and 
minimum  curves  correspond  to  values  obtained  from  the  curves  of  Fig. 

3-2  at  an  altitude  of  125  n.m. 

For  this  typical  case,  with  an  effective  target  velocity  of  25,230 
feet  per  second,  an  uncertainty  in  critical  frequency  of  1  Me  will  pre¬ 
vent  accuracies  of  better  than  1  foot  per  second  at  operating  frequencies 
below  about  350  Me  and  greater  than  0.1  foot  per  second  below  about 
1100  Me.  Increasing  the  critical  frequency  accuracies  by  a  factor  of 
10  increases  the  velocity  accuracy  by  the  same  factor. 


Figure  6-7  is  an  elaboration  of  the  information  in  Fig.  6-6  for 
the  particular  case  of  f^  ^  400  Me.  Doppler  velocity  error  is  shown 
as  a  function  of  uncertainty  in  critical  frequency  and  the  maximum, 
mean,  and  minimum  curves  are  shown. 
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